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Formation of stationary localized states in one- 
dimensional chain due to a linear impurity and a modi- 
fied nonlinear impurity is studied. Furthermore, a one- 
dimensional chain with linear and modified nonlinear site 
energies at the alternate sites is studied and rich phase dia- 
grams of SL states are obtained for all the systems we con- 
sidered. The results are compared with those of the linear 
and nonlinear systems. 



I. INTRODUCTION 

It is well known that the transport properties of a 
system is directly related to the formation of stationary 
localized (SL) states in the system. Localized states ap- 
pear due to the presence of impurity or disorder (which 
breaks the translational symmetry) in the system H. 
There has been detailed study on the formation of SL 
states in various systems in the presence of linear impu- 
rities. One the other hand the effect on the formation of 
localized states due to the presence of nonlinear impuri- 
ties has been studied recently by the use of the discrete 
nonlinear Schrodinger equation given as [|]-[l6) 



dC n 
i—— 
dt 



e n C n + V(C n+1 + Cn-O + X n\C n \ a C n , (1) 



where a is equal to 2 or arbitrary constant, C n is the 
probability amplitude of the particle to be at the site 
n, e„ and \n are the static site energy and the nonlin- 
ear strength at site n respectively, and V is the nearest 
neighbor hopping element. The nonlinear term {Cn^Cn 
arises due to interaction of the exciton with the lattice 
vibration [Jl7],[l8|. If the lattice oscillators are purely 
harmonic in nature, then a = 2. The Eq. (1) has been 
derived with the assumption that the lattice oscillators 
are much faster in motion compared to that of the exci- 
ton and the lattice oscillators are independent of others. 
Equation (1) has been used to study the formation of 
SL states due to a single nonlinear impurity and that 
due to a dimeric nonlinear impurity in one dimensional 
chain as well as in a Cayley tree. A perfect nonlinear 
chain with nonlinear impurity have also been considered 
to study the SL states and it has been shown that SL 
states appear even thou gh t he translational symmetry of 
the system is preserved |K]|-|lq| . Very recently the effect 
of the presence of both the linear and nonlinear impuri- 
ties in one dimensional chain and in the Cayley tree has 
been studied and rich phase diagrams of SL states have 
been obtained The Eq. (1) may be modified if the 
lattice oscillators are assumed to be coupled with their 



nearest neighbors. Thus, under the assumption that the 
oscillators are coupled to its nearest neighbors and are 
much faster in motion compared to the motion of the 
exciton, the discrete nonlinear Schrodinger equation is 
modified and takes the form: 



dC n 
i—— 
dt 



£nC n + V(C n +i + C n -l) 



+ Xn (\C n+1 \ 2 + |C„_!| 2 + 2\C n \ 2 )C n 



(2) 



where C n , e n , V and \n carries same meaning as for the 
Eq. (1). The modified discrete nonlinear Schrodinger 
equation (Eq. (2)) has been derived from the coupled 
exciton-oscillator system by Kopidakis et al. I^Ofl . We 
notice that the Eq. (2) contains more nonlinear terms 
compared to Eq. (1) and, in addition, Eq. (2) is also rel- 
evant in condensed matter physics because the nearest 
neighbor coupling among the lattice oscillators is taken 
into account in deriving the equation. One interesting 
property of Eq. (2) is the self-trapping property and this 
aspect has been studied by Kolosakas et al. |2l|. An- 
other interesting property of this equation is that it can 
produce SL states. In fact the Eq. (2) has been exploited 
to study the formation of SL states in one dimensional 
system in the presence of a single nonlinear impurity, 
and that in the presence of a dimeric nonlinear impu- 
rity [p2| . The Eq. (2) has also been used to find the SL 
states in a perfectly nonlinear chain j22|. However the 
system with both the linear and the modified nonlinear 
impurities has not been considered yet to study the for- 
mation of SL states. Thus we are interested to observe 
the mixed role of linear and modified nonlinear impu- 
rities on the formation of SL states in one dimensional 
system. 

The organization of the paper is as follows. The one 
dimensional chain with a linear impurity and a modified 
nonlinear impurity is considered in Sec. I. In Sec. II 
we consider the one dimensional chain with linear and 
modified nonlinear impurities at alternate sites of the 
lattice. Finally In Sec. Ill we summarize our findings. 



II. ONE DIMENSIONAL CHAIN WITH A 
LINEAR IMPURITY AND A MODIFIED 
NONLINEAR IMPURITY 

Consider a one-dimensional chain consisting of a lin- 
ear impurity of strength e placed at the zeroth site and 
a nonlinear impurity of strength x at the first site. The 
relevant Hamiltonian for such a system can be written 
as 



c*c n 



+1, 



n— — oo 



1 



-f (\C2\ 



|C | 2 + 2|C 1 | 2 ) ICJ 2 , 



(3) 



where the nearest-neighbor hopping matrix element is 
assumed to be unity. As a possible solution for the sta- 
tionary localized states, we consider 



iBi 



(4) 



where E is the stationary state energy. It is also well 
known that impurity states in one-dimensional system 
are exponentially localized. Therefore, the presence of 
two impurities at the two consecutive sites (site and 
site 1) in the chain, suggests us to consider the following 
forms of <b„: 



n-l 



n > 1 



K = <i>\ [sgn(£')r/] 
0_ W = <Msgn(£)77] H n<0. 



(5) 



The form of 4> n given above are exact and may be de- 
rived from the Green's function analysis M . Here 77 lies 
between and 1 and is given as 



1 

V= 2 



\E\ - VE 2 -A 



(6) 



Since we are dealing with the localized states, (f>o and <j>\ 
can be assumed to be real without forsaking the math- 
ematical rigor. Three possible cases may be considered, 
i.e., 4>x = 4>q, 4>i — —0o, and <fii ^ <f>o, the solutions for 
which are called, respectively, symmetric, antisymmet- 
ric, and asymmetric solutions. To take into account of 
all three possibilities, we introduce a variable (3 = 4>o /4>i ■ 
Consequently, (3 = 1 corresponds to the symmetric so- 
lutions, (3 = —1 to the antisymmetric solutions, and 
/3 7^ ±1 to the asymmetric solutions. Substituting the 
ansatz of Eq. (5) and the definition of (3 to the normal- 
ization condition ^2 \C n \ 2 = lj we obtain 



\<h 



1 



if 



1 + /3 2 



(7) 



Using Eqs. (4), (5) and (7) in Eq. (3), we obtain the 
effective Hamiltonian for the reduced dynamical system: 



H cfi = sgn(E)r) H 



(1 + /? 2 ) [ ' 



4 (l + /3 2 )< 



(/3V + 2/3 4 + /? 2 ). 



(8) 



We first consider the case of (3 = ±1. The effective 
Hamiltonian for this case becomes 



H± = sgn(E) V ± (L_^) + 1(1 - r, 2 ) 
+ A(3 + r 7 2 )(l-r; 2 ) 2 , 



(9) 



where "+" refers to the symmetric states and "— " cor- 
responds to the anti-symmetric states. For a particular 
system, e and \ are constants, while 77 is a variable which 
determines the energy of the localized state. The value 
of 77 should be determined self-consistently and, for this 



reason, we treat 77 as the dynamic variable in the ef- 
fective Hamiltonian of the reduced dynamical system. 
The fixed-point solutions corresponding to the localized 
states can be found, from the variational principle, by 
the condition dH^ s /dr) — with 77 S [0,1]. Since the 
form of the probability profile is rigorous W, we obtain 
the correct values of 77, satisfying the equation 



8 _ 77(1 — 77 2 )(3i7 2 + 5) 
^± ~ sgn(£) - (| ± 1)7? 



f±(v)- 



(10) 



From Eq. (10) it is clear that the states will appear 
above the band for x > and below the band for 
X* 1 < 0. [Note that the band edges are at E/V = ±2.] 
We will consider positive values of x an d, hence, the 
sgn(E) may be taken to be positive. For symmetric 
states, the value of f+(rf) exhibits the extreme values: 
f+{rj) — * as 77 — > and f+(r]) — * 00 as 77 — > l/(e/2 + l). 
Thus, f+(r)) increases as 77 increases from 77 = and 
diverges at 77 = l/(e/2 + 1), and it becomes negative 
for?7 > l/(e/2+l). This means that there will al- 
ways be only one solution (symmetric state) for all pos- 
itive values of e and \- For antisymmetric states, we 
found that f-(rj) — > as 77 — > and f-{rj) — * 00 as 
77 — » l/(e/2 — 1). The divergence of f-(rf) is, however, 
confined between and 1 if and only if e > 4. There- 
fore, for e > 4, f-(rj) increases as 77 increases from 77 = 0, 
diverges at 77 = l/(e/2 — 1), and finally becomes nega- 
tive for 77 > l/(e/2 — 1). This implies that there will 
always be one antisymmetric state for e > 4 for any 
value of x > 0. However, the situation is different for 
e < 4; the value of /-(??) approaches as 77 — » and 
77 — > 1 and remains positive and finite for < 77 < 1. 
This implies that f-{rj) has at least one maximum in 
the range of < 77 < 1 . We confirm by a graphical anal- 
ysis that there is only one maximum. Therefore, there 
will be a critical value of x separating two antisymmet- 
ric states from the region of no antisymmetric state for 
e < 4. The dashed line in Fig. 1 is such a critical line 
and is denoted by Xcr- Thus, there are three regions in 
the entire plane of e versus Xi having different number 
of symmetric and antisymmetric SL states. The region 
of e > 4 has one symmetric and one antisymmetric SL 
states, while the region of e < 4 is divided into two re- 
gions by the critical line x%-> below of which has only 
one symmetric state and above of which has one sym- 
metric state and two antisymmetric states. Therefore, 
the maximum number of SL states contributed by the 
symmetric and anti-symmetric solutions is three. 

Next, we consider the asymmetric solutions, i.e., for 
(3 7^ ±1. The effective Hamiltonian in Eq. (8) should be 
considered as a function of two variables, 77 and (3, while 
e and x remain constant for a particular system. The 
fixed-point solutions corresponding to the SL states can 
be obtained from the following two extreme conditions 



dH eff /d v = 0, 
dH ee /d[3 = 0. 



(11) 



After a trite algebra, the fixed-point equations reduce 
to 



2 



2(1 + /3 2 )(l + /3 2 - 277/3 -e?7) 

-xn(i - v 2 ){WW + 3/3 4 + 2d 2 ) = 

and 

2(l + /3 2 )(l-/3 2 -e/3) 
+X0Q- -V 2 ){2(3 2 v 2 + 3/3 2 + 1) = (12) 

These coupled equations can not be solved analytically 
and, accordingly, we analyze them numerically to find 
the number of asymmetric SL states (number of solu- 
tions with r\ G [0, 1]). 
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FIG. 1. The phase diagram for SL states in the plane of 
e versus x f° r a system of a one-dimensional chain with a 
linear impurity and a modified nonlinear impurity. Different 
regions are marked with the number of states and the kind 
of states; the region marked by "ls+2a+2as" means that 
the region has one symmetric, two antisymmetric and two 
asymmetric SL states. The small unmarked region has only 
one symmetric SL state. 

It is found that there are three critical lines sepa- 
rating different number of asymmetric SL states in the 
plane of e versus x denoted by the solid curves in Fig. 1. 
Maximum number of asymmetric SL states obtained is 
two. Considering all possible states we find that max- 
imum number of SL states is five. The "s", "a" and 
"as" (in Fig. 1) stand for symmetric, antisymmetric 
and asymmetric states, respectively. For example, the 
region marked as "ls+2a-|-2as" has one symmetric SL 
state, two antisymmetric SL states, and two asymmetric 
SL states. It should be noted that the number of states 
changes one by one as we go from one region to another 
in the phase space. As an example, as we cross from the 
region "ls+las" to the region "ls+2a+las" , the num- 
ber of states increases by two; however, the number of 
states on the critical line which separates two regions is 
"ls+la+las" , yielding the increment of the number of 
states one by one. The maximum number of SL states 
in one dimensional chain with two linear impurities was 
found to be two Q and that in one dimensional chain 
with two modified nonlinear impurities was found to be 
four . We thus conclude that one linear impurity and 
one modified nonlinear impurity play a role to increase 
the maximum number of SL states. The phase diagrams 



of SL states here is very complicated when compared to 
the that for one dimensional chain with two modified 
nonlinear impurities. 

Figure 2 shows the variation of energy of the states as 
a function of x f° r a fixed value of e = 6. We see that 
three SL states appear at x = 5 as expected from the 
phase diagram (Fig. 1). The energy of the symmetric 
state (top solid line in Fig. 2) and that of the anti- 
symmetric state (another solid line in Fig. 2) increases 
while the energy of the asymmetric state decreases as 
X increases. At x ~ 6.06 one more asymmetric state 
appear whose energy increases with the increment of 
X and finally both the asymmetric states disappear at 
X — 6.55. Thus we find that the symmetric state and the 
anti symmetric state are strongly localized compared to 
the asymmetric states in the parameter regime consid- 
ered in Fig. 2. Through stability analysis [p2| it can be 
shown that the state whose energy decreases with the 
increment of x is unstable. 
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FIG. 2. Variation of energy as a function of \ IS shown for 
e = 6. 



III. LINEAR AND MODIFIED NONLINEAR 
SITE ENERGIES AT ALTERNATE SITES 

In this section, we consider a one-dimensional chain 
with alternate site energies; the even sites have a linear 
energy e and the odd sites have a modified nonlinear 
energy x{\C n +i\ 2 + |C„_i| 2 + 2|C„| 2 ). Such a system 
can be described by the Hamiltonian given as 

1 v x 

H= 2 (^n^n+l + C*C n+ i) + - |C 2n | 2 

n— — oo n— — oo 

oo 

+ f £ (|C 2n +2| 2 + |C 2 „| 2 + 2|C 2 „ +1 | 2 )|C 2li+1 | 2 . (13) 

n— — oo 

For x — 0, the system reduces to a periodically modu- 
lated linear system. The translational symmetry of the 
system is preserved for both x — and x / 0. Us- 
ing C n = <j) n exp(-iEt), we obtain the Hamiltonian in 
terms of <f> n . It is, however, not possible to find the exact 
ansatz for the localized states in this system, although a 



3 



certain rational ansatz can be considered. For example, 
the on-site peaked solution, inter-site peaked or inter- 
site dipped solutions are possible. We will consider each 
of these in the subsequent subsections. 
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FIG. 3. The phase diagram of SL states in the e versus 
X plane for the system of a one dimensional chain with lin- 
ear and modified nonlinear site energies at alternate sites. 
The number '2' added in every region indicates that there 
are two on-site peaked states apart from the symmetric, an- 
tisymmetric and asymmetric states; one of them is peaked 
at the linear site and another one is peaked at the nonlinear 
site. For example, the region marked by 'ls+2a+las-|-2' has 
one symmetric state, two antisymmetric states, one asym- 
metric state and two on-site peaked states (one peaked at 
the linear site and another one peaked at the nonlinear site). 



A. On-site peaked solution 

Since the system has linear and nonlinear sites at reg- 
ular intervals, the solution may have a peak either at 
the linear site or at the nonlinear site. We first consider 
the solution peaked at the linear site, i.e. the solution 
is peaked at the central site of the system. We therefore 



consider the ansatz 



M 



Using this and the 



normalization condition, we obtain the effective Hamil- 
tonian as 



H e ff = 



2 V 



(l + r? 4 ) 
I + 77 2 ' 2 (1 + 77 2 ) 2 



+ 



x v 2 (i-v 2 ) 

2 (1 + 77 2 )(1 + r7 4 ) ' 



(14) 



The fixed-point equation, dH c g/dr) = 0, reduces to 

2(1 - 77 8 )(1 + ?7 4 ) - 2677(1 - /7 2 )(1 + V 4 ) 2 
+ X r ] ((l + 7 1 2 ) 2 (l-2r ] 2 )(l + r ] i ) 



-77 2 (l-r7 4 )(l + 277 2 + 377 4 )) = 0. 



(15) 



Solving Eq. (15), we found that there is always one SL 
state peaked at the linear site for any value of e and X- 
The localized solution peaked at the nonlinear site may 
be obtained considering that each site is shifted by one 
lattice site along either positive or negative direction, 
and the resulting fixed-point equation is given as 



2 (1 + r/ 2 )(l - V 4 )(l - r/ 8 ) 2 + 2o?(l - r/ 4 )(l - t? 8 ) 2 

+ XV ((1 + V 2 )(l V 8 f(l 3r, 2 ) + 8r, 6 (l - ry 4 ) 2 ) 

- XV (4(1 - r/ 4 )(l - V le ) + 3r?(l - r, 2 )(l - V s ) 2 ) = 0. (16) 

It is found from the analysis of Eq. (16) that only one SL 
state peaked at the nonlinear site appears irrespective 
of the values of \ an d £■ 



B. inter-site peaked and dipped solutions 

For the inter-site peaked and dipped solutions we use 
the dimeric ansatz in Eq. (5). After short calculations, 
we obtain the effective Hamiltonian for the reduced dy- 
namical system as 



H eff = sgn(E)r] + 



+ 



P(l-V 2 ) e (1+/3V) 
1 + /3 2 2 (1 + /3 2 )(1+t? 2 ) 

x(i-v 2 ) 2 h 2 (i + /? 4 ) 



4 (1 + /3- 



2\2 



(1 



2^2 



+ 



x (i - v 2 ) 

4(1 + /? 2 ) 2 



V°) 
(P 4 



+ /3 2 

+ v 4 ) 



r? 6 (l + /3 4 ) 
(1 — ?7 8 ) x (l-rr 8 ) 



(17) 



with (5 and 77 as defined earlier. We first consider the 
symmetric and antisymmetric states, i.e., (3 = ±1. The 
fixed-point equation yields 



4 _ r7((l + ^ 2 )(l + 2r? 2 -v 4 ) -V 4 ) 



X 3 



(lTr?)(l + r/ 2 ) 2 



= f ± (v), (18) 



where "+" refers to the symmetric states and "-" to 
the antisymmetric states. It should be noted that e has 
no role in the formation of symmetric and antisymmet- 
ric SL states in the system. From Eq. (18) it is clear 
that only one symmetric state appears for any finite 
value of X- However, for antisymmetric states, the sit- 
uation is different. The value of f~(v) in Eq. (18) ini- 
tially increases as 77 increases and finally decreases as 
77 approaches to 1, leaving a maximum at r\ G [0,1]. 
Thus, there are two critical values of x sucn that there 
is no anti-symmetric SL state for X < 9.07, two anti- 
symmetric SL states for 9.07 < X < 10.66 and one anti- 
symmetric SL state for X > 10.66, shown as the two 
vertical lines in Fig. 3. 

We next consider the asymmetric states of j3 7^ 1. In 
this case e plays its role, unlike the case of symmetric 
and antisymmetric states. The fixed-point equations are 
given as 

(1 - t? 8 ) 2 [(1 + ry 2 ) 2 (l + /3 2 )(1 + [3 2 - 2(3 V ) - er?(l - /3 4 )] 

+xv(±-v 4 ) 2 (i-v & ) 

1(1 + /3 4 )(1 + 3,? 4 ) - /3 2 (1 + ry 2 )(l + ry 4 ) + rf 



+X^ 3 (l-^)(l- ? 7 8 )(l + /3 4 ) 
^-(l + ^Xl + ry 4 )] 

-xv(i - v 4 ?(P 4 + v 4 ){i + v 2 + v 4 -v 6 ) = o 



(19) 



and 



4 



(l-^Xi-^Xi-/? 4 ) 
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- X/ 3(l-r 7 2 )(l-ry 4 )[r 7 2 (l + /3 4 )(l + ?? 4 ) 
+/3 4 (l-77 8 )+/3 4 + r7 4 ] =0 



(20) 



Solving Eqs. (19) and (20) we obtain two critical curves 
(solid curves in Fig. 3) separating the \ — e plane 
into various regions with various number of asymmet- 
ric SL states as shown in Fig. 3. Symmetric states are 
denoted by 's', antisymmetric states are denoted by 'a' 
and asymmetric states are denoted by 'as'. The region 
marked by 'ls+2a+las+2' contains one symmetric SL 
state, two antisymmetric SL states, one asymmetric SL 
state, one SL state peaked on the nonlinear site and 
one SL state peaked at the linear site. The maximum 
number of states obtained in this case is 6 which is the 
same as for the chain with modified nonlinear impuri- 
ties of different strengths at alternate sites. Thus in this 
case the linear impurity plays the same role as that of 
nonlinear impurity to obtain the maximum number of 
SL states. No SL state is obtained for x = 0. This is 
expected because the system with \ = reduces to a 
periodically modulated linear system which again may 
be renormalized to a perfect linear system. 



IV. SUMMARY 

Formation of stationary localized states in one dimen- 
sional chain due to the presence of a linear impurity 
and a modified nonlinear impurity is studied. Maximum 
number of SL states in this case is found to be five which 
is larger than that obtained for a one dimensional chain 
with two modified nonlinear impurities. Thus one may 
conclude that the presence of both linear and nonlinear 
impurities plays a role to increase the maximum number 
of SL states. The variation of energy of SL states as a 
function of the nonlinear strength for a fixed value of 
linear impurity strength is presented. It is found that 
the symmetric and antisymmetric states are strongly lo- 
calized when compared with the asymmetric states. 

Furthermore, a one dimensional chain with linear and 
modified nonlinear energies at alternate sites is consid- 
ered. Here the maximum number of SL states is found 
to be six. The linear energy strength at alternate sites in 
the chain does not play any role to produce the symmet- 
ric and antisymmetric SL states and the system does not 
produce any SL state for x=0 which is expected. How- 
ever, surprisingly, the linear energy strength plays a role 
to produce the asymmetric SL states. Rich phase dia- 
grams of SL states for both the systems considered here 
are presented. 
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